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Abstract 

We calculate the cross sections for dissociation of J / ip by pions and p- mesons 
within the framework of a meson exchange model. We find that these cross 
sections are small at center-of-mass energies less than 1 GeV above threshold, 
and that dissociation rates are less than 0.01 c/fm in a thermal meson gas at 
temperatures where such a description makes sense. 
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A pronounced suppression of the production of strongly bound heavy quarkonium states 
(J/ip, T) is considered as a signature for the formation of a new, deconfmed phase of 
strongly interacting matter in collisions of heavy nuclei at high energy, the quark-gluon 
plasma |L|. Experimental investigations of charmonium production in nuclear reactions 
have been carried out for over a decade at the CERN-SPS (p+A, O+U, S+U, Pb+Pb) 
and the Fermilab Tevatron (p+A) These studies have shown that the J/ip and ip' 

production on nuclear targets is, indeed suppressed relative to expectations from nucleon- 
nucleon reactions. 

For p+A collisions the observed suppression of J / ip and ip' can be explained as absorption 
of a common precursor, probably a non-resonant, color-octet (cc) state, on nucleons ||. 
This systematics extends to O+U and S+U reactions for the J/if>, but not for the ip', 
where additional suppression is observed ||. This effect can be quantitatively described 
as absorption of the ip' on comoving secondary hadrons ( "comovers" ) , mostly light mesons 
]7] [1IJ. Recently, measurements of charmonium production in the Pb+Pb system [[J have 
revealed the presence of an additional "anomalous" suppression mechanism also for the J/ip. 
It is presently under debate whether this mechanism can be absorption of the J/ip or the 
down-feeding Xc on hadronic comovers, or whether a more drastic explanation, such as quark 
deconffnement, is required 111-113 . 



In spite of the obvious need to understand the physics of charmonium absorption on 
comoving hadrons, only very few quantitative predictions for J/ip absorption cross sections 
on light mesons have been published [|14],|15|]., and these predictions differ by several orders 



of magnitude. We hope that the results presented below will help clarifying this situation. 
Our calculation is based on an effective hadronic Lagrangian describing the interactions 
among 7T, p, J/ip, D, and D* and, for the sake of simplicity, on the assumption of thermal 
equilibrium among the light hadrons. 

The most abundant mesons in a hot hadronic gas are ir, K, and p. In addition, nucleons 
may be abundant at high baryon chemical potential. In the following, we will not consider 
kaon induced reactions, but the framework presented here can be easily extended to include 



this absorption channel. At low energy, 7r and p mesons, as well as nucleons, can induce 
dissociation processes when encountering a J/ip particle (we list only the reactions with the 
lowest energy thresholds): 

7T + J/V> -> D + D*, D + D*; (1) 
p + J /if, -> D + D; (2) 
N + J/ip A c + D. (3) 

The kinematic thresholds for these reactions are (|l|) 640 MeV, (fj) —135 MeV, and (H) 115 
MeV, respectively. The reaction p + J/ip — > D + D has the lowest total invariant mass 
threshold and is, in fact, exothermic. At thermal equilibrium it is only suppressed due to 
the relatively high mass of the p-meson which causes p-mesons to be less abundant than 
pions in a thermal hadron gas. The reaction (||) is mostly of interest in the fragmentation 
region, where the nucleon density is high. We will not consider it further here, because it is 
expected to be of less importance for J / ip suppression at central rapidity 

There are two possible approaches to the J / ip dissociation problem: one at the quark level 
where the large mass of the c-quark is used to separate perturbative from nonperturbative 
aspects of the problem, and another one that makes use of effective hadronic interactions. 



The former approach is based on the pioneering work of Peskin |L6| and Bhanot and Peskin 



1 1711 , w ho realized that interactions between heavy quark bound states and light hadrons 
can be described perturbatively, if the heavy quark mass tuq is sufficiently large. The small 
size of the (QQ) state allows for a systematic multipole expansion of its interaction with 
external glue fields, where the color-dipole interaction dominates at long range. 

In the Bhanot-Peskin framework, light hadrons interact with the J/ip only via their 
glue content. Kharzeev and Satz |TJ| applied this formalism to inelastic ir — J/ip scattering 
and showed that the absorptive cross section is proportional to the component of the gluon 
structure function Gh{x) of the interacting hadron which is sufficiently energetic to dissociate 
the J/ip. Because of this requirement of rather hard gluons, only highly energetic hadrons 
are capable of exciting a J/ip above the dissociation threshold. The pion momentum must 



reach 5 GeV/c before attaining an absorption cross section in excess of 1 mb. The thermally 
averaged cross section, in this framework, remains less than 0.1 mb for a pion gas within 
any realistic temperature range. 

One limitation of the Bhanot-Peskin approach lies in the fact that the J/ip is not truly 
a Coulombic bound state but probes also the confining part of the cc-potential. Although 
these nonperturbative contributions to J /if} absorption by light mesons can be estimated 
18| , the reliability of such estimates is difficult to assess. 

From a microscopic point of view, the reactions (|l|J^) can be viewed as quark exchanges, 
where the J / if) transmits a charm quark to the light meson and picks up a light (u, d, or s) 
quark. Since similar reactions among light hadrons typically have large cross sections at mod- 
erate energies, one may suspect that these reactions also proceed with significant strength 
above their respective kinematic thresholds. It is therefore of interest to calculate the disso- 
ciation cross sections in the framework of an essentially nonperturbative approach based on 
a hadronic model that incorporates quark confinement. The cross section for the reaction 
(0) was calculated in the nonrelativistic quark model by Martins, Blaschke, and Quack |T5 



in the first Born approximation. Including D*D, DD*, and D*D* final states, the total 
cross section was found to peak around 1 GeV energy above threshold (in the center-of-mass 
system) at a value of about 7 mb. However, the large magnitude of their result is due to 
the action of the long-ranged, confining interaction between the quarks, which is modeled as 
a "color-blind" attractive interaction between the quarks with a Gaussian momentum de- 
pendence. Because this interaction is taken as attractive independent of the color quantum 
numbers of the affected quark pair, it does not cancel in the interaction of a light quark 
with a point-like (cc) pair in a color-singlet state, in contrast to the one-gluon exchange 
interaction. 

Since the results obtained within the constituent quark model differ by orders of magni- 
tude from those obtained in the Bhanot-Peskin approach, and because they depend critically 
on the particular implementation of the confining quark-quark interaction, it seems prudent 
to calculate the charm-exchange cross sections in a second, entirely different framework of 



hadronic reactions. In the effective meson theory, the exchange of a (cq) or (cq) pair, where 
q, q stands for any light quark, can be described as the exchange of a D or D* meson between 
the J/ip and the incident light meson. 

At the hadronic level, the charm exchange reaction between the J/tp and a light hadron 
can proceed either by exchange of a D- or a /}*-meson. In fact, Regge theory dictates that 
the charm exchange reaction is dominated by the exchange of the _D*-trajectory in the high 
energy limit. However, here we are not interested in charm exchange at high energies but near 
the kinematical threshold, because the relative motion of the hadrons is limited to thermal 
momenta. We shall only consider the D-exchange reactions. (The D*-exchange reaction 
cross sections exhibit an unphysical rise with energy due to the exchange of longitudinally 
polarized D*-mesons, which would need to be eliminated with the help of form factors.) 

In order to calculate the various Feynman diagrams for the reactions we need to 

construct the effective three-meson vertices. We can do so by invoking a strongly broken 
local U(4) flavor symmetry with the vector mesons playing the role of quasi-gauge bosons. 
Denoting the 16-plet of pseudoscalar mesons (n, i], r]', K, D, D s , rj c ) by $ = <piTi, where the Tj 
are the U (4) generators, and the vector meson 16-plet (u, p, <fi, K*, D*, D*, t])) by V M = V^Ti, 
the free meson Lagrangian reads 

Co = tr(^$ f <9 M $) - tr^V^X^V, - 9„V M ) - tr^MpS) + -tr(V^M v V M ). (4) 

Here Mp and My denote the mass matrices for the pseudoscalar and vector mesons, respec- 
tively. Because of the heavy mass of the charm quark, Mp and My break the U(4) symmetry 
strongly down to U(3), the mass of the strange quark introduces a weaker breaking to U(2), 
and the axial anomaly further breaks the U(2) symmetry to SU(2) in the case of the pseu- 
doscalar mesons. All these symmetry breakings are embodied in the physical mass matrices. 
It is convenient, in the following, to work with the mass eigenstates. 

The meson couplings are obtained by replacing the derivatives <9 M by the "gauge covari- 
ant" derivatives = <9 M — igV^. In first order in the coupling constant g this procedure 
leads to the following interactions 
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£ int = ig tr (Vy^$ _ + ig tr (d^ v [V^ V v \ - [V^, V^K) . (5) 

If the U(4) flavor symmetry were exact, we would expect all couplings given by the same 
constant g. In view of the significant breaking of the flavor symmetry we anticipate that the 
effective coupling constants for different 3-meson vertices will have different values. However, 
we will see below that the coupling have remarkably similar values. 

In order to describe it- and p- induced J/tp dissociation, we need the ipDD, ttDD*, and 
pDD vertices. From (g) we derive the following interactions: 

£nDD* = T^gnDD* (pTiD^d^i - DTiD*^ - h.C.) 

£ p dd = ^9 P DDPi [pTid^D - d^DnD} 

jC^ dd = ig^DDV (DdpD - (d„D)D) . (6) 

Here the Tj denote the generators of SU(2)-isospin symmetry in the fundamental represen- 
tation (Pauli matrices). The coupling constants g^DD, and g P DD can be derived from the 
.D-meson electric form factor in the standard framework of the vector meson dominance 



model [19]. If 7v denotes the photon- vector meson V mixing amplitude, the usual analysis 



yields the relations 

2 

IpQpDD = em 2 p , I^dd = ^ em l ( 7 ) 

The photon mixing amplitudes 7v can be determined from the leptonic vector meson decay 
widths: r Vee = a7v/(3rav). Inserting the experimental numbers, one finds 

g P DD ~ 5.6, g-ipDD ~ 7.7. (8) 

As is well known, different ways of deriving the value of the p-meson coupling yield values 
differing by about 20%. The coupling constants (§) must therefore be considered to be given 
with an error of this order of magnitude. 

The 7r-meson coupling between D and D* can be obtained from the decay width of the 
Z}*-meson: D* — > Dn. Unfortunately, only an upper bound for this decay rate is known at 
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present [pO| , corresponding to g n DD* < 15. A theoretical estimate for this decay rate, based 



on QCD sum rules [plf| , gives 

9-irDD* « 8.8. (9) 

For lack of an experimentally determined value of g n DD* we adopt this value here. 

Analytical expressions for the amplitudes and cross sections for the processes correspond- 
ing to the Feynman diagrams in Figure [I] are easily derived. We begin with the two diagrams 
(|I|a,b) for absorption by pions. Averaging over initial and summing over final isospin and 
spin states, one obtains 

2 2 \ 2 \ / ( I m 2 *yri 2 \ 2 



ml 



-zHdd9-kDd*7^ 4m x 2 ^ m D 

6 w t 2 \ m 2 D * ) \ 



(10) 



Here s, t, and u are the Mandelstam variables and we have introduced the notation t' = 
t — m 2 D) u' = u — m 2 D . The contributions of the diagrams (|l]c,d) need to be added coherently, 
because they lead to identical final states. However, their results are related by crossing 
symmetry (t *->■ u). After some algebra one finds: 



\M c + M d \ 2 = ^gl DD g 2 DD 



l(w D -^)(w D -M^; 



+— Us - A,nl - H-»H-^)\ ( 2s _ Aml _ 
u't' \ mi \ m A 



+ (i! «-> u') 



(11) 



The charm exchange cross sections are plotted in Figure ^| as functions of the center-of- 
mass energy y/s. The tt + J/ip cross section (dashed line) starts at zero, because the reaction 
is endothermic, whereas the p+J/tp cross section (solid line) is finite at the threshold because 
of its exothermic nature. Note that the cross sections for these two J/tp absorption reactions 
are of similar magnitude over the energy range relevant to a thermal meson environment. 

Although pions are far more abundant than p-mesons in the temperature range where a 
hadronic gas is likely to exist, most of these pions do not have sufficient energy to dissociate 



a J/if). If one counts only those pions having energy above the kinematical threshold for 
J/"0-dissociation, n- and p-mesons are about equally rare, because their effective density not 
exceeding 0.1 fm -3 even at T = 200 MeV, where the hadron picture of a thermal environment 
probably already breaks down. As a consequence, we expect dissociation rates of similar size 
for 7r and p mesons at thermal equilibrium. Figure |3| shows these J/ ip absorption rates as a 
function of temperature. Even at the (unrealistically) high temperature T = 300 MeV, the 
thermal dissociation rates are still so small that a J /if} continues to have a lifetime around 
10 fm/c. For T < 200 MeV, dissociation by hadrons is completely negligible on the time 
scale of the lifetime of a hot hadronic gas state in nuclear collisions. 

The inclusion of a form factor in the meson vertices of the Feynman diagrams ([I]) would 
further reduce the dissociation rates. E.g. a Gaussian form factor exp(— Q 2 /Ql) with Q = 
1.5 GeV would reduce the absorption rates by slightly more than one order of magnitude. 
This implies that our result represents an upper limit for the hadronic dissociation rate. 
As we noted above, our analysis remains incomplete because of the neglect of D*-exchange 
reactions. It would be most interesting to include these reactions in the framework of a 
complete effective Lagrangian describing the interactions of the pseudoscalar and vector 
mesons with quark content (qq), (Qq), (qQ), and (QQ), where q stands for any light quark 
flavor and Q denotes a heavy quark. Such a Lagrangian embodying chiral symmetry for the 



light quarks has recently been given by Chan | 22| . Loop corrections to the tree diagrams 
would permit the study of form factor effects in this approach, as well. 

In conclusion, we have calculated the dissociation cross sections and thermal absorption 
rates of J/ if> mesons on pions and p-mesons in the framework of an effective meson exchange 
model. We find that the dissociation rates in a thermal meson gas, at temperatures where 
such a gas is expected to exist, are very small corresponding to a survival time of the J/if> in 
excess of 100 fm/c. Our results are in line with those recently obtained for other absorption 



channels, such as n + J/if> — > ir + if) 1 and p + J/ip^-K + rjc ||24|| . Taken together, these 
results indicate that a thermal hadronic environment in the confined phase of QCD cannot 
be the cause of significant absorption of J/if> mesons. 
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FIG. 1. Lowest order Feynman diagrams contributing to the charm exchange reactions (||) and 
(|2|). Single lines represent pseudoscalar mesons; double lines denote vector mesons. 
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FIG. 2. Cross sections for the charm exchange reactions described by the diagrams of Figure 
|], as functions of cm. energy. Dashed line: pions, solid line: />mesons. 
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FIG. 3. Thermal J/tp absorption rates as function of temperature. The pion and p-meson rates 
are shown separately. 
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